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Abstract 



We investigate in the equal-time formalism the derivation and truncation 
of infinite hierarchies of equations of motion for the energy moments of the 
covariant Wigner function. From these hierarchies we then extract kinetic 
equations for the physical distribution functions which are related to low- 
order energy moments, and show how to determine the higher order moments 
in terms of these lowest order ones. We apply the general formalism to scalar 
and spinor QED with classical background fields and compare with the results 

derived from the three-dimensional Wigner transformation method. 
PACS: 03.65.Bz, 05.60,+w, 52.60.+h. 
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I. INTRODUCTION 



Transport theory [|TJ based on the Wigner operator is extensively used to describe the 
formation and evolution of highly excited nuclear matter produced in relativistic heavy ion 
collisions. The Wigner operator can be defined in 4-dimensional or 3-dimensional []7]|| 
momentum space, which we denote by VV(x,p) and W(x, p), respectively Correspondingly, 
there are two formulations for the phase-space structure of any field. Either of these two 
formulations has its advantages and disadvantages. Besides its manifest Lorentz covariance 
(which is very useful from a technical point of view), another characteristic feature of the 
4-dimensional formulation for QCD @-|5| and QED || is that the quadratic kinetic equation 
can be split up naturally into a transport and a constraint equation. The complementarity of 
these two ingredients is essential for a physical understanding of quantum kinetic theory |§ . 
In the classical limit, these two equations reduce to the Vlasov and mass-shell equations, 
respectively. The main advantage of the 3-dimensional approach f7||| is that it is easier 
to set up as an initial value problem: one can directly compute the initial value of the 
Wigner operator from the corresponding field operators at the same time. In the covariant 
frame this is not possible since the covariant Wigner operator is defined as a 4-dimensional 
Wigner transform of the density matrix and thus includes an integration over time. Hence 
in this approach the initial condition for the Wigner operator at very early times must be 
constructed phenomenologically. Some true quantum problems like pair production fH| in 
a strong external field have thus so far been solved only in the 3-dimensional (or equal-time) 
formulation 0JT1]. 

One way |7|] to obtain equal-time kinetic equations which parallels the procedure in the 
covariant formulation is to Wigner transform the equation of motion for the equal-time 
density operator $(cc, y, t). For spinor QED this procedure results in the BGR equations |7j 
for the equal-time Wigner functions. In Ref. Jl2|| we suggested a different derivation which 
is based on taking the energy average of the covariant kinetic equations in the 4-dimensional 
formulation. It exploits the fact that the equal-time Wigner function is the energy average 



2 



(i.e. zeroth order energy moment) of the covariant one. With this method we showed 
for spinor QED that the direct energy average of the covariant kinetic equations leads, in 
addition to the BGR transport equations for the spinor components of the equal-time Wigner 
function, also to a second group of constraint equations which couple the equal-time Wigner 
function to the first order energy moment of the covariant one. In the classical (H — > 0) limit, 
these additional equations provide essential constraints on the equal-time Wigner function 
and allow to reduce the number of independent distribution functions by a factor of two 



12| . In the general quantum case, the additional equations determine the time evolution 
of the energy distribution function which in general can not be expressed in terms of the 
equal-time Wigner function. In this sense the BGR equations do not provide a complete set 
of equal-time kinetic equations. 

As we will discuss in this paper, this incompleteness has a more general aspect. As just 



mentioned, the equal-time Wigner operator is related to the covariant one by [12 



W(x,p) = I dEW(x,p,E) , (1.1) 

where we wrote p = (E,p), E independent of p. As such it is only the lowest member of an 
infinite hierarchy of energy moments of the covariant Wigner operator: 

W J (x,p) = JdEE j W(x,p,E), j = 0,l,2,---, (1.2) 

with W (x,p) = W(x,p). Therefore, to set up a complete equal-time transport theory 
which contains the same amount of information as the covariant theory one needs dynamical 
equations for all the energy moments. Any covariant kinetic equation will thus correspond 
to an infinite hierarchy of coupled kinetic equations for its energy moments, i.e. for the 
equal-time Wigner operators Wj(x,p). 

This infinite hierarchy only exists for genuine quantum problems where the energy can 
exhibit quantum fluctuations. In the classical limit, the covariant Wigner operator satisfies 
the mass-shell constraint p 2 = E 2 — p 2 = m 2 , and the energy dependence of the covariant 



Wigner function thus degenerates to two delta- functions at E — ±E P = ±^/m 2 + p 2 . The 



equal-time Wigner operator W(x,p) then splits into a positive and a negative frequency 
component, 

W(x,p) = W + (x,p) + W~(x,p) , (ft->0) (1.3) 

and all energy moments can be expressed algebraically in terms of these as 

W J (x,p) = W;(x,p) + W-(x,p) 

= E{,W + (x,p) + (-E p yW-(x,p), j = 0,l,2,... (fc-O). (1.4) 

The solution of the equal-time kinetic equations for W(x,p) thus also determines the dy- 
namics of all higher energy moments. Thus, in the classical limit, a simple zeroth order 
energy average of the covariant kinetic theory yields a complete equal-time kinetic theory. 

In the general quantum case, the higher order energy moments Wj(x,p), j > 1, contain 
genuine additional information and can no longer be expressed algebraically through the 
equal-time Wigner operator W(x,p). This means that in principle in the equal-time formu- 
lation we are stuck with the problem of solving an infinite hierarchy of coupled equations. 
Actually, there are two such hierarchies, one resulting from the covariant transport equa- 
tion ("transport hierarchy"), the other arising from the generalized mass-shell constraint 
("constraint hierarchy"). In practice this raises the problem of truncating the hierarchy in 
a physically sensible way. Since only the low-order energy moments of the covariant Wigner 
function have an intuitive physical interpretation, it turns out that physics itself suggests an 
appropriate truncation scheme. We will show that the hierarchies of moment equations are 
structured in such a way that the first few low-order moments form a finite and closed sub- 
group of equations which can be solved as an initial value problem, and that (surprisingly) 
all the higher order moments can be derived from these low-order moments recursively using 
only the constraint hierarchy, i.e. without solving any additional equations of motion. The 
equations from the transport hierarchy for the higher order moments are redundant. 

We first discuss on a general basis, starting from the covariant approach, the deriva- 
tion and truncation of equal-time hierarchies of kinetic equations. For illustration we then 
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consider in full generality the case of a transport theory for scalar fields with arbitrary 
scalar potentials. For this case everything can be worked out explicitly to arbitrary order 
of the moments. We give the subgroup of equations which fully characterize the first few 
low-order moments, prove the independence and redundance of the transport equations for 
the higher order moments outside this subgroup, and obtain from the constraint hierarchy 
explicit expressions for all the higher order moments in terms of the solutions of the low- 
order subgroup. We then apply the general formalism to scalar and spinor QED. Here the 
equations have a more complicated structure, and we restrict our attention to the closed 
subgroup of equations for the lowest order moments, discussing the redundance of the trans- 
port equations for the higher order moments and their recursive determination through the 
constraint hierarchy only for the first moments outside the closed subgroup of equations 
for the low-order moments. We will compare our results with the kinetic equations for the 
equal-time distribution functions obtained previously in Refs. [f7|,p|,[12| . Our final result will 
be a complete set of kinetic equations which can be implemented numerically as an initial 
value problem. 



II. GENERAL FORMALISM 

The 4-dimensional Wigner transform of the equation of motion for the covariant density 
operator leads to a complex, Lorentz covariant kinetic equation for the Wigner operator. It 
couples the one-body Wigner operator to two-body correlations Q, which in turn satisfy 
an equation which couples them to three-body terms, and so on. After taking an ensemble 
average this generates the so-called BBGKY hierarchy |13[ for the n-body Wigner func- 



tions. A popular way to get a closed kinetic equation for the one-body Wigner function 
(i.e. the ensemble average of the one-body Wigner operator) is to truncate the BBGKY hi- 
erarchy at the one-body level, by factorizing the two-body Wigner functions in the Hartree 
approximation. So far most applications of quantum transport theory have employed this 
approximation, and in the following we will also restrict ourselves to it. For us the mean 
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field approximation provides a crucial simplification, and at present it is not obvious to us 
how to generalize our results in order to include correlations and collision terms. 

For a scalar field in mean field approximation the complex equation for the self-adjoint 



scalar Wigner function can be separated into two independent real equations fl2f : 

G(x,p)W(x,p) = 0; (2.1a) 
F(x,p)W{x,p) = 0. (2.1b) 

The first equation corresponds to a generalized Vlasov equation; after performing the en- 
ergy average it generates a hierarchy of transport equations for the energy moments Wj (x, p) 
("transport hierarchy"). The second equation is a generalized mass-shell constraint; it gen- 
erates a hierarchy of non-dynamic constraint equations ( "constraint hierarchy" ) . Equations 
with the structure given in ( |2.1|) will be the starting point for our discussion of scalar field 
theories in Sees. |TXT| and |IV A| . Factors of p M in the dynamical operators G(x,p) and F(x,p) 



arise from the Wigner transformation of the partial derivative in the Klein-Gordon equa- 
tion. Since the latter contains at most second order time derivatives, at most two powers 
of po occur. In fact, G(x,p) is linear in p while F(x,p) is quadratic in p . This will be 
important below (see Sec. [11 B| ). 

For spinor fields the covariant Wigner function is a 4 x 4 matrix in spinor space which 
for a physical interpretation must be decomposed into its 16 spinor components W s . In this 
way the complex kinetic equation for the Wigner function matrix is split into 32 independent 
real equations for the self-adjoint spinor components |6|JT2|. These equations can be further 
divided into two subgroups according to their anticipated structure after performing the 
energy average: 



Hi 



G ss ' (x, p) W s ' (x, p) = , (2.2a) 

s'=l 
16 

£ F ss '(x,p) W s '(x,p) = , (s = 1, 2, . . . , 16). (2.2b) 



s'=l 



The first subgroup leads to equations containing only first order time derivatives and thus 
generates 16 hierarchies of transport equations for the energy moments of the spinor com- 
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ponents ("transport hierarchies"). The other subgroup which involves both first and second 
order time derivatives leads to a set of 16 hierarchies of constraint equations ("constraint 
hierarchies") for the equal-time moments of the spinor components. For the lowest energy 
moments, the spinor components W s (x, p) of the equal-time Wigner function, the details 
of this procedure were worked out in Ref. [O, and we will use these results in Sec. [IV C 



Since the original Dirac equation is linear in the time derivative, the dynamical operators 
G ss ' (x,p) and F ss '(x,p) contain at most single powers of po- in fact, the operators G ss '(x,p) 
are independent of po- 



A. Hierarchy of energy moments 

In this subsection we will concentrate for simplicity on a single covariant kinetic equation 
of the generic form 

G(x,p)W(x,p) = 0, (2.3) 

where G(x,p) contains at most two powers of po and of the space-time derivative operator 
<9 M , but an arbitrary number of derivatives with respect to the momentum space coordinates 
(see Appendix^. We will return to the full set of equations (|2.1|) resp. (|2.2|) in the following 
Sections. 

We begin by decomposing the energy dependence of the Wigner function into a basis of 
orthogonal polynomials hj(E): 

oo 

W(x,p,E) = ^2w j (x,p)h j (E) . (2.4) 

The expansion coefficients Wj(x,p) are defined in the equal-time phase space. Using the 
orthonormality relation 

J d f i(E)h i (E)h j (E) = 5 ij , (2.5) 

where dfi(E) is the appropriate integration measure associated with the chosen set of poly- 
nomials hj, the equal-time components Wj(x,p) can be related to energy moments of the 
covariant Wigner function constructed with the basis functions hj(E): 
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Wj(x,p) = J dn(E)hj(E)W(x,p,E) . (2.6) 

If the system has finite total energy the covariant spinor components must vanish in the 
limit E — > ±00. We will assume that they vanish at infinite energy faster than any power 
of E such that for any combination of integers i,j,m,n > 0, we have 

I ME) dE {^( h + E ) Em )^-m^P,E))=0. (2.7) 

With exponential accuracy we may therefore restrict the energy integration to a finite interval 
—A < E < A. Introducing the scaled energy u = E/A we can thus use as our set of basis 
functions the Legendre polynomials 



h n {u) = ^±ip„( w ) (2.8) 

with the trivial measure d/j,(cj) = du on the interval [—1,1]. 

As discussed above, the dynamical operator G(x,p) in ( |2.3p in general contains powers 
of E up to second order and an infinite number of energy derivatives d/dE. In terms of the 
new dimensionless energy variable we may thus write 

M 00 / o \ n 

G(x,p,u)= £ Y.G mn (x lP )u m — , (2.9) 

m=0n=0 \ aUJ J 

with M < 2. Substituting this double expansion into equation ( |2.3| ), multiplying by hi(ou) 
from the left and integrating over energy uj we obtain 

£G mn (x,p)/ dwhi{u)u m — ;W(x,p,u) = 0. (2.10) 

m,n J - 1 

Inserting the expansion ( |2.4j) of the covariant Wigner function W(x, p, u), this can be written 



as 



where 



J2H lJ (x,p)w J (x,p) = 0, i = 0,1,2,..., (2.11) 

3=0 



M 00 

Hv(x, P ) = J2E C T d rnn(x,p) (2.12) 
m=0 n=0 
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with 



C™ n = f 1 duoh^^dlh^uo) . (2.13) 



It is easy to see from ( |2.13| ) that 



C™ = for n>j and i > j + m - n . (2.14) 

For j > i — m + n the coefficients are in general nonzero. Therefore, the sum over j in 
Eq. (12111) extends over the range j > max[0,i — M\. For each value of i, Eq. ( |2.11| ) thus 



contains an infinite number of terms. In this form Eqs. ( p.ll| ) are thus not practically 
useful. However, one can use the surface condition ( [2.7|) to rewrite Eqs. (|2.11| ) in such a 
way that each equation contains only a finite number of terms. Returning to Eq. ( |2.10| ) and 
integrating by parts, we can replace the integrand by 

h l (uj)uj m dZW(x,p 1 uj)= (2.15) 

n-l 

E("R W^V ffl )0(^,w)) + (-) n ^(h i Hw m )w(i, P)W ) . 
z=o 

The contribution to the integral in fl2.10| ) from each term in the sum is fully cancelled by 
the surface condition (|2.7|) , and only the last term in (|2.15|) survives. Inserting it into ([2.1 0|) 
and using again the expansion ( p.4| ) we find instead of Eqs. (|2.11| - pT3"D the following set of 
equations: 

i+M 

9ij(x,P)wj(x,P) = 0, i = 0,1,2,..., (2.16) 

3=0 



with 



M i+m 

g^p) = E E c^G mn (x,p) (2.17) 

m=0 n=0 



and the coefficients 



= i > /(2t + l)(2j + l) £ i duP j (cj)(-&r(P i (u,)u, m ) . (2.18) 

The latter can be determined recursively from = 5^ (which results from the orthogonality 
relation ( |2.5j )) by using the recursion relations for the Legendre polynomials, see Appendix 
Since the nonvanishing coefficients c™ n are now restricted to the domain 

9 



n < i + m and j < i + m — n , (2-19) 

the sum over j in (|2.16|) runs now only over the finite range < j < i + M. The first 
inequality in ( |2.19| ) was already used in ( |2.17| ) to limit the sum over n. 

The Pi{uS) are polynomials in uj of order i, and thus the equal-time components Wj(x,p) 
occurring in Eqs. fl2.16p are linear combinations of the energy moments Wk(x,p) of order 
k < i (see Eq. (|1.2|) ). For each value of i, Eq. (|2.16|) thus provides a relation among the 
first i + M + 1 energy moments of the covariant Wigner function W(x,p,uj) (including 
the zeroth order moment W(x,p) = v2wo(x,p)). As i is allowed to run over all positive 
integers, Eqs. ( [2.16|) form an infinite hierarchy of relations among the energy moments of the 



covariant Wigner function. Each covariant equation of the type (|2.3|) generates its own such 
hierarchy. Only the full set of these infinite hierarchies of moment equations constitutes a 
complete equal-time kinetic description of the system under study. 

B. Truncating the hierarchy 

In order to discuss possible truncation schemes we must return to the complete set of 
covariant kinetic equations. Let us concentrate here on the scalar case, Eqs. (|2.1|) , and write 
down the two resulting hierarchies of moment equations as 

i+M 

9ij{x,p)w j {x,p) =0, (2.20a) 

3=0 
i+M+1 

Y fepK(i,p) = I (i = 0,1,2,...). (2.20b) 

3=0 



In writing down the upper limits of the sums we already used that F(x,p) in ( p. lap contains 



one power of po more than G(x,p) in ( |2.1a ). For the scalar field case one has M = 1. For 



the spinor case one obtains from ( |2.2| ) a similar set of equations with M = 0. 

Let us now try to truncate these hierarchies for the moments Wj at some order j max . 
The equations from the "transport hierarchy" (|2.20a|) with hierarchy index % < I t involve all 



moments Wj with < j < I t + M, i.e. the first I t + M + 1 moments (including the lowest 
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moment with index 0). Similarly, the equations from the "constraint hierarchy" ( |2.20b|) with 



hierarchy index i < I c involve the first I c + M + 2 moments < j < I c + M + 1. For a 
closed set of equations both hierarchies must be truncated at the same order j max , i.e. we 
must have 

I t + M = I c + M+ l= j max . (2.21) 

Truncating in this way we are left with Jj+1 equations from the transport hierarchy and J c +1 
equations from the constraint hierarchy. In order solve them the number of equations must 
at least equal the number of moments. However, if there are more equations than moments, 
the system may be overdetermined, and therefore we would like to require equality of the 
number of equations and moments: 

I t + 1 + I c + 1 = j max + 1 . (2.22) 

The two conditions ( |2.2ip and ( 2.22 ) have a unique solution: 



I t = I c + l = M, Jmax = 2M, (2.23) 

which yields M + 1 transport and M constraint equations for the first 2M + 1 energy 
moments. Smaller values of j max don't yield enough equations, and larger values lead to 
an (at least superficially) overdetermined system of equations. For spinor fields (M = 0) 
the truncated set involves only one transport and no constraint equation; for each spinor 
component it gives a single kinetic equation (the BGR equation 0) for its lowest energy 
moment, the equal-time Wigner function W s (x,p). For scalar fields (M = 1) the truncated 
set contains two transport equations and one constraint for the three lowest order moments 
w , Wi, w 2 . 

If we go beyond this minimal closed subset of moments and equations, we get two equa- 
tions for every additional moment, one from the transport hierarchy and one from the 
constraint hierarchy. As we will show explicitly in the next Section, in the constraint hier- 
archy (|2.20b| ) the highest moment always comes with a constant coefficient. As we increase 
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the hierarchy index % in Eqs. (|2.20|) , at each step the newly occurring moment can thus be 
explicitly expressed in terms of the already known lower order moments using the corre- 
sponding constraint equation from ( |2.20b| ) . As we will discuss, these higher order constraint 
equations contain important physics. But in addition, at each step there is also a dynamical 
equation of motion for the new moment from the transport hierarchy ( |2.20a| ). How can the 
two equations be consistent? The answer is that this transport equation is not an indepen- 
dent new equation, but (with some algebraic effort) can be expressed as a combination of the 
lower order equations which have already been used. Our proof of this fact uses explicitly 
the structure of the dynamical operators G mn and F mn . It involves cumbersome algebra, 
and only for scalar fields with only scalar potential or mean field interactions we have been 
able to find a general proof. For scalar and spinor QED the proof is still incomplete, and 
we will only demonstrate the first step for the 2M + 2 nd moment. A completion of the proof 
presumably requires a so far missing deeper insight into the general dynamic structure of 
the moment equations and their relation to the underlying covariant theory. 



III. SCALAR FIELD THEORY 

In this Section we will give an explicit and complete discussion of the moment hierarchy 
for the simplest case of a scalar field theory in Hartree approximation. We exemplify the 
truncation of the hierarchy and the recursive computation of the higher order moments 
beyond minimal truncation. The discussion in the following Section for the practically more 
relevant case of QED will be technically more involved and, unfortunately, also less complete. 

A. Covariant kinetic equations 

Consider the Klein-Gordon equation with a scalar potential U(x): 

(dl + ml + U{x)) 4>{x) = . (3.1) 
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The covariant Wigner function is the four-dimensional Wigner transform of the covariant 
density matrix y) = (<&(x,y)): 

W(ar, p) = J d 4 y e ip ' v ${x, y) = J d 4 y e ip ' v (x + §) 0+ (x - §) ) . (3.2) 

To derive the kinetic equations for the scalar Wigner function, we calculate the second-order 
derivatives Q<9^ + d^j and (|<9^ — d*j\ of the covariant density operator, and then employ 
the Klein-Gordon equation ( O ) and its adjoint. After taking the ensemble average and 
performing the Wigner transform we obtain two complex kinetic equations: 

(R 2 - p 2 + m 2 + U (x - id p ) - ip-d x ) W{x, p) = , (3.3a) 
(\dl -p 2 + m 2 + u(x + id p ) + ip-8 x ) W{x,p) = . (3.3b) 

Since the scalar Wigner function is real, adding and subtracting these two complex equations 
yields two real equations of the type ( J2.1[ ) . After reinstating % the corresponding operators 
G and F are given by 

G(x,p) = hp-d x + lmM 2 (x,p) , (3.4a) 

h 2 

F(x,p) = -p 2 + — d 2 x + ReM 2 (x,p) , (3.4b) 
where the mass operator M 2 is defined as 

M 2 (x,p) =m 2 + T, e (x,p) + i£, (x,p) , (3.5a) 

t, e (x,p) = cos (~Y^J U ( x ) > (3.5b) 

t Q (x,p) = sin f-^-J u ( x ) ■ ( 3 - 5c ) 

Here the triangle operator A is defined as A = d x -d p where the coordinate derivative d x 
acts only on the scalar potential U(x) while d p acts only on the Wigner function. 



B. Semiclassical expansion 



In the general quantum situation the particles have no definite mass due to quantum 
fluctuations around their classical mass shell and collision effects in the medium. In the 
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situation here with only an external potential this is illustrated by the mass operator M 2 . 
Only in the classical limit % — > it reduces to the quasiparticle mass 

ReM 2 (i,p) = m 2 (x) = m 2 + U{x) , (3.6a) 
lmM 2 {x,p)=0. (3.6b) 

In this case the constraint equation reduces to the on-shell condition 

(p 2 - m 2 {x)) W (x,p) = (3.7) 

for the classical covariant Wigner function Wo- The classical transport equation arises from 
the general transport equation at first order in %. The first order contribution to the mass 
operator is 

ReM 2 (x,p) = 0, 

1m M 2 (x,p) = hm(x) (d x m(x))-d p , (3-8) 

and we obtain the covariant Vlasov equation 

(p-d x + m(x)(d x m(x))-dp)W (x,p) = , (3.9) 

with a Vlasov force term induced by x-dependent effective mass term. For scalar fields there 
is no first order quantum correction to the operator F in ( |3.4| ), and from the zeroth order 
term we obtain the mass-shell condition for the first order Wigner function: 

(p 2 -m 2 (x)) Wi{x,p) = . (3.10) 

This discussion holds universally for arbitrary potentials. If, for instance, U(x) is generated 
by the scalar field (j)(x) itself in Hartree approximation, 

U( X ) = - C + \(4>(x)4> + (x)} = -C + X J ^YV(x : p), (3.11) 

with a mass parameter C and a coupling strength A, this model provides a useful tool for a 
dynamical description of spontaneous symmetry breaking |§. 

14 



C. The 3-dimensional dynamical operators 



We now perform the energy average of the covariant transport and constraint equations 
( j2.ip and construct the hierarchy (|2.20| ) of moment equations. The first step is the double 
expansion of the type ( |2.9|) for the covariant dynamical operators G(x,p) and F(x,p): 



Here 



Gmn(Xi Pi 



1, n 







hAdt for m 

hp-'Vx — o for m = n = 

— (^r^o) for m = 0, n ^ even 

-^Gtr^e) for m = 0,nodd 

else 



{x,p,t) 



-A 2 



for m = 2, n = 



^-d, 2 ; + p 2 + m 2 , + a e for m 



n 







MI) 

-^(f) n (^o) 



for m = 0,n^0 even 
for m = 0, n odd 
else. 



(3.12a) 



(3.12b) 



a e (x,p,t) = cos ^-V^-VpJ f/(a;,t) 
a (a;,p,t) = sin f^V x -V p J U(x,t) 



(3.13a) 
(3.13b) 



are the three-dimensional analogues of the covariant operators S e and S in Eq. ( |3.5| ). 
Again, the spatial gradients act only on U(x,t), while the momentum gradients act on the 
equal-time Wigner functions (i.e. on the energy moments Wj(x,p)). 

The three-dimensional dynamical operators G mn (x,p) and F mn (x, p) must now be com- 
bined with the coefficients c™ n to obtain the dynamical operators §ij(x,p) and fij(x,p) 
which are needed in the transport and constraint hierarchies. This is done in Appendix [B]. 
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D. Minimal truncation 



The resulting transport hierarchy is truncated at It = M = 1, the constraint hierarchy 
at I c — M — 1 = 0. This yields the following equations for w (x,p), Wi(x,p), and w 2 (x, p): 



goowo + goiwi = , (3.14a) 
9ww + gnw 1 + g 12 w 2 = , (3.14b) 
/oo^o + /oi^i + fo2W 2 = . (3.14c) 

The dynamical operators and are given in Appendix and Eqs. (|3.12p . Reexpressing 
Wj in terms of the energy moments Wj from Eq. ( |1.2| ), 



w (x, p) = -^=W(x, p, t) , (3.15a) 
Wl {x,p) = ^^|wi(x,p,t) , (3.15b) 
Mx,P) = ^Jl(^W 2 (x,p,t)-W(x,p,tj) , (3.15c) 



Eqs. Q3.14 ) can be rewritten as 



d t Wi{x,p,t) = - (p-V x - ~a (x,p,t)\ W(x,p,t) , (3.16a) 
d t W 2 (x, p,t) = - (p-V x - ~&o(x, P, tj\ Wt(x, P,t) + ~ (d t a e {x, p, t)) W(x, p, t) , (3.16b) 
W 2 (x,p,t) = ^ (d 2 t - Vj) +p 2 + m 2 Q + & e (x,p,t?j W(x,p,t) . (3.16c) 

Please note that all powers of the cutoff A cancel in the final expressions as they should. 

The two transport equations ( |3.16| a,b) do not decouple, not even in the classical limit 
h — *■ 0. To achieve decoupling one must return to the covariant equations in Sec. |111 A| and 



study their semiclassical limit as given in Sec. [llllj| before performing the energy average. 
Then the mass-shell condition ( |3.7| ) can be used to rewrite all higher order energy moments 
in terms of the zeroth order moment as explained in the Introduction, Eq. ( |1.4j ). With this 
information the constraint (|3.16cj ), in the limit h — ► 0, becomes trivial, 
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W 2 (x, p, t) = E 2 (x, t) W(x, p, t) , E 2 (x, t)=p 2 + m 2 (x, t) , (3.17) 
while the two transport equations ( [3.16a| ) and (|3.16b| ) become identical and can be written 



in the form of a Vlasov equation for the charge density (see Sec. |V|): 



dtW x {x,p,t) + ( ~ V x E p ■ V p j W 1 (x iP ,t) =0. (3.18) 

The reason why the information contained in Eq. ( |1.4j ) cannot be easily recovered directly 
from the 3-dimensional transport and constraint equations is that in their derivation, through 
Eq. (p.!5| ), we made heavy use of partial integration with respect to the energy. In the 
classical limit this has the unfortunate effect of spreading the information contained in the 
on-shell condition over the whole infinite hierarchy of 3-dimensional constraint equations. 

Although we have always talked about Eq. ( |3. 16c|) as a "constraint equation" it is clear 
that, as far as solving the minimal subset ( |3.16|) of equal-time kinetic equations is concerned, 
this terminology is only adequate in the classical limit h — > 0. In general quantum situations 
it is a second order partial differential equation for the lowest order moment W(x, p, t) which 
must be solved as an initial value problem together with the first order partial differential 
equations (|3.16a|) and (|3.16b|). 



E. Higher order moment equations 

The next higher moment u>3 is determined by the third equation in the transport hier- 
archy and the second equation in the constraint hierarchy, 

^20^0 + #21^1 + ^22^2 + <?23^3 = , (3.19a) 

fww + fuwt + fi 2 w 2 + /13W3 = . (3.19b) 
With the dynamical operators from Appendix and Eqs. (|3.12|) and 

w 3 (x,p) = ±^(JLw 3 (x,p,t) - ~Wi(x,p,t)) (3.20) 
we obtain (using Eq. ( |3.16a| )) 
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d t W 3 (x,p,t) = - (p-V x - -a (x,p,t) s j W 2 {x,p,t) + (d t a e (x,p,t)) W x {x,p,t) 



--(d 2 t a {x,p,t))W{x,p,t) 



W 3 (x, p, t) = fa (d 2 t - V 2 ) +p 2 + m 2 + a e (x, p, t)j W 1 (x, p, t) 

h 

-- (d t a (x,p,t))W(x,p,t) . 



(3.21a) 



(3.21b) 



By substituting (|3.19b| ) into (|3.19a|) and taking into account the following commutators: 



[Gio, Fqo] — — 2G01-F20 j 



[G 



00, r oo\ 



02-^20 ) 

K 2 h 2 h 2 

[—d 2 , a e/o ] = —dcr e/o -d + —d 2 a. 



r e/o ) 



[p-V x , a e/o ] = p-V x a e/o + -V x a o/e -V x 

h 2 

[P 2 , v e /o\ = -hp-V x a o/e + —V 2 a e/o , 



(3.22a) 
(3.22b) 
(3.22c) 
(3.22d) 

(3.22e) 



the third transport equation ( |3.19a| ) can be rewritten in terms of the first transport equation 
( glial) as 



foo (#00^0 +0OlWl) = 0. 



(3.23) 



This implies that the transport equation ( |3.19a|) for ^3 is redundant. The third-order mo- 



ment W3 is completely determined in terms of the solutions of the minimal subgroup (|3.16|) 
by the constraint equation (|3.21b| ). It arises from Eq. ( |3.19b|) by noting that /13 is a constant, 



/13 = C3 = — |^|A 2 (see Appendix |B|) , and solving for u> 3 : 



w 3 



~T (flQWo + fllWi + f 12 w 2 ) ■ 



(3.24) 



Note that, in contrast to (|3.16cj ), Eq. (|3.21b|) does not require solving a partial differential 
equation because everything on the r.h.s. is known from the solution of Eqs. ( |3.16| ). 

The above procedure can be extended to all higher order moment equations, by repeat- 
edly using the commutators listed in ( |3.22| ). In general one finds that a transport equation 



E 9ij W j = 
3=0 



(3.25) 



with i > 2 can be re-expressed in terms of the first (i — 1) transport equations as 

i-2 /j+l \ 

£/*-iU \Y,9ikWk) =0. (3.26) 

3=0 \k=0 J 

Thus, except for the first two, all transport equations are redundant. The higher order 
moments W{ with i > 2 can be computed from their constraint equations 

1 - 

w i = -77H^-2,j w j, (3-27) 
°* i=o 

with the constant C{ given by 

Ci = fi-2,i = f ' = • (3.28) 

(2« — 1)^(2* — 3)(2i + 1) 



IV. APPLICATION TO QED 
A. Scalar QED 

In this Section we discuss the application of the general formalism developed in Sec. |I|to 
QED. Since some of the equations will be rather lengthy we will economize on the notation 
by dropping all factors of H. The latter are correctly given in Refs. JT2]JT7[] to which we refer 
in case of need. 



We begin with the case of scalar QED with external electromagnetic fields. In Ref. [12 
we discussed the semiclassical transport equations for this theory by energy averaging the 
semiclassical limit of the covariant transport equations. In this subsection we will derive the 
general equal-time quantum transport equations by performing the energy average without 
any approximations. In the following subsection the result will be compared with the cor- 
responding equations derived by directly Wigner-transforming the equations of motion for 
the equal-time density matrix ||. 

In scalar QED the scalar field obeys the Klein-Gordon equation 

+ ieA^(x))(d^ + ieA»(x)) + m 2 ) </>(x) = , (4.1) 
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Following an analogous procedure as in Sec. [Ill A| (see [12| for details) one derives two 
covariant kinetic equations of type ( |2.1| ) for the covariant Wigner function 



W(x,p) = ( dV iM (|(x + f) 



exp 



1/2 

i< I ds A(x + sy)-y 

1/2 



0+ (x - |) 



The corresponding covariant dynamical operators G and F are given by 
G(x,p) = fl"(x,p)t)^x,p) , 

F(x,p) = ^D f "(x,p)ID fl (x,p) - Uf i (x,p)U. li (x,p) + m 2 , 



(4.2) 

(4.3a) 
(4.3b) 

tl^(x,p) =p il -ie l^dss F^(x - isd p ) d u p , (4.3c) 

D^(x,p) =d^-e / * dsF^ v (x - isd p ) d v v . (4.3d) 

■piiv _ qiij^v _ Q"A^ i s the electromagnetic field tensor. 

The structure of the equal-time transport theory for scalar QED is very similar to that 
for a scalar potential U(x) which we considered in the previous section. The difference 
resides solely in the expressions for the dynamical operators G mn and F mn . In Appendix 
we provide the double expansions (|2.9|) for the basic operators IP and -D M as well as for G 
and F. 

In terms of the covariant Wigner function the charge current j M and the energy- 
momentum tensor T^ v are expressed as 



f[x,p) = ep^ W(x,p) , 
T^(x,p) =p^p u \V(x,p). 



(4.4a) 
(4.4b) 



After performing the energy average these equations translate into relations between the 
first three moments Wq,Wi,W2 and the equal-time phase-space distributions for the scalar 
density W(x,p, t), the charge density p(x,p,t), and the energy density e(x,p,t): 

1 



w (x,p,t) 



wi(x,p,t) 



ks/2 

fs 1 



W(x,p,t) 



2eA 2 



p(x,p,t) 



w 2 (x,p,t) = 2^^/~ (j^e(x,p,t) - W(x,p,t)^j . 



(4.5a) 
(4.5b) 
(4.5c) 
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The charge current density j(x,p, t) and the momentum density P(x,p, t) can be expressed 
in terms of W and p as 

P 

j(x,p,t) = epW(x,p,t), P(x,p,t) = -p(x,p,t) . (4.6) 



The subgroup ( |3.14| ) determining the moments wq, w\ and W2 can thus be equivalently 



rewritten as transport and constraint equations for W, p and e: 

irfiP + (d t 7r + Tz-d)W = , (4.7a) 

d t e + -(d t 7t + ird + D)p + {tt D + tt-J - d t A + G-d)W = , (4.7b) 

e 

e = (±{d\ - dt) - (tt 2 - tt 2 ) + m 2 + A) W - 2 -^p , (4.7c) 

The expressions of g,ij and were obtained via the relations given in Appendix ^ from the 
equal-time operators ( |C2| ) in Appendix y. We have also used the commutators 

[dt, ttq] = D , (4.8a) 
[4 i] = . (4.8b) 

Due to the line integrals over s in the operators 11^ and which guarantee the gauge 
invariance of the formalism the equal-time operators G mn and F mn for QED are much 
more complicated than Eqs. ( [3.12|) for the case of a scalar potential. This is the origin of the 



more complicated structure of the minimal subgroup (|4.7| ) of 3-dimensional kinetic equations. 
Please note the sequence of the operators in Eqs. (|4.7| ): in particular the generalized time 
derivatives dt act on everything following them. Eqs. ( |4.7|) are thus much more intricately 
coupled than the corresponding equations ( |3.16| ) for the pure scalar case. 

It is instructive to integrate Eqs. ( |4.7|) over p to obtain equations of motion for the 
corresponding space-time densities: 

d t p(x,t) + V-j(x,t) = 0, (4.9a) 
d t e(x, t) + V-P{x, t) - E(x, t)-j(x, t) = , (4.9b) 
e(x,t) = J^^(d?-V^+p* + m^W(x,p,t). (4.9c) 
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The first two equations express the conservation of electric charge and of energy-momentum 
while the last equation gives the energy density in terms of the scalar equal-time 
phase-space density W(x,p,t) including quantum corrections. In Eq. ( |4.9b[ ) P(x,t) = 
f(d 3 p/(2n) 3 )P(x,p,t) is the momentum density in coordinate space, and the last term de- 
scribes the conversion of field energy into mechanical energy by the work done by the electric 
field on moving charges. 

Let us now consider the next two equations in the hierarchy, the transport and constraint 
equations ( [3.19|) for the third-order moment W3. Again the constraint equation can be 



directly solved for W3 in terms of the solutions wq,wi, w 2 of Eqs. (|4.5| ), (|4.7|) (see Eq. (|3.24|) ). 
Our task is to show that the transport equation for u> 3 is redundant. To this end we substitute 
the constraint ( |3.19b|) into the transport equation ( |3.19a| ) and use the scalar QED analogue 



of the commutators ( |3.22|) , namely ( |4.8| ) and 



[4 B) = 3F, (4.10a) 

[4 tt] = -I+V x 7t Q , (4.10b) 

[J, tt] =-V x -G, (4.10c) 

[d, 7T ] = V X 7T , (4.10d) 

[d t , d] = -eV x -V p (\ ds E[x + isV p , t) , (4.10e) 

[tt, 7r ] = -ie j^dss E(x + isV p , t) + eV x / ^ ds s 2 E(x + isV p , t)-V p , (4.10f) 

as well as the identity 

J\ds (2is + ^V X -V P - s 2 V x -V p ^j E(x + isV p ,t) =0 , (4.11) 

to rewrite ( |3. 19a| ) in the form 

foo (Qoo^o + Qoiwi) + foi {9wWq + g n wx + gi 2 w 2 ) = . (4.12) 

Note that in the derivation of the commutators (|4.10|) we used Maxwell's equation 

= , (4.13) 
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for the dual field strength tensor = ^e^ uljp F ap . The identity ( |4.11| ) (which has no 
analogue in a theory without gauge invariance) is proven by expanding the electric field 
E(x + isV p ,t) = e ls ^ x ' T7p E{x ) t) and integrating term by term. Eq. ( |4.12|) expresses the 
transport equation for w% in terms of the transport equations for W\ and Wi- This proves 
that it is redundant. 

^From the comparison of (|4.12| ) with (|3.23| ) we observe that in the case of scalar QED the 



third transport equation is expressed in terms of the first and second transport equations, 
while only the first one occurs in the case of scalar potentials. This difference results from the 
line integrals which occur in the gauge theory. In a transport theory with gauge invariance, 
the kinetic momentum is not but fl M . The energy average of the zeroth component of the 
second term in ( [4.3c|) yields the equal-time operator tc given in (|C2|) , and n = — Ir/oi i n 



turn generates the coefficient in front of the second bracket in Q4.12 ). In a transport theory 



without gauge freedom all coefficients /j^+i vanish due to the absence of linear terms in p M 
in the covariant constraint equation. 

We have not had the patience to carry the above considerations to higher orders in the 
energy moments. The corresponding calculations become extremely messy. Based on the 
experience with pure scalar theory we expect all higher order transport equations to be 
redundant, but we have failed to discover a simple calculational technique which permits us 
to prove this in an elegant way. 

B. Comparison with the Feshbach-Villars approach 

In Ref. H a set of equal-time transport equations for scalar QED was derived by directly 
Wigner transforming the equations of motion for the (non-covariant) equal-time density 
matrix. Since the Klein-Gordon equation ( |4.1| ) contains a second-order time derivative, 
its direct translation into 3-dimensional phase space does not lead to a sensible transport 
equation which should have only first-order time derivatives. Therefore the 3-dimensional 
approach exploits the Feshbach-Villars representation |14| of the field equations of motion 
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which contains only first-order time derivatives. The price to pay (in addition to the loss 
of manifest Lorentz covariance) is the introduction of an auxiliary field which results in a 
rather complicated 2x2 matrix structure of the scalar Wigner function. For the energy 
averaging method advocated here there is no such problem. In this case second-order time 
derivatives appear only in the constraint hierarchy, and the transport hierarchy contains 
only first-order time derivatives. 

Expressing the equal-time density operator <&(x,y,t) of the scalar field in terms of the 
Feshbach-Villars spinors, performing a 3-dimensional Wigner transformation and using the 
Klein-Gordon equation in Feshbach-Villars form, the authors of || obtained the follow- 
ing coupled equations for the 4 spinor components fi(x,p,t) (i = 0, 1,2,3) of the Wigner 
function: 

mdtfo = -4t-d(f 2 + / 3 ) , (4.14a) 
md t h = -(id 2 - tt 2 )(/ 2 + / 3 ) + 2m 2 / 2 , (4.14b) 
mdth = (^d 2 - tt 2 )/i + 7r-df - 2m 2 h , (4.14c) 
md t f 3 = -(^d 2 - 7z 2 )U - 7i-d/ . (4.14d) 

In Ref. || these equal-time phase-space distributions were related to the physical phase-space 
densities p, j, e, and P as follows [3, |l5|1 : 

P = e/ , (4.15a) 

e=^(/ 2 + /3)+m/ 3 , (4.15b) 

3 = %(h + fs) , (4.15c) 

P = pf . (4.15d) 



Unlike the spinor decomposition for spinor QED |6|,|7] where each component of the Wigner 
function has a definite physical meaning, there is no obvious physical interpretation for the 
component f\ in scalar QED. 

The comparison of the charge currents j in Eqs. (|4.15c ) and ( O ) leads directly to the 
identification 
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W=-(f 2 + h). (4.16) 



Using the relations ( |4.15| ) and (|4.16| ) one finds from a suitable linear combination of the 



kinetic equations (|4.14j ) the following equations of motion for the physical phase-space dis- 
tributions p and e, 

-d t p + Tt-dW = , (4.17a) 

2 



d t e + ~*-d,p - ~ f - tt 2 )^ + rfVJ W = , (4.17b) 
e = - ^d 2 ) + i(7T 2 + p 2 ) + m 2 ) W , (4.17c) 



2 I v 4 

J 2 \ 1 

<*) + -. 

.4^ 1 2 > 2 

as well as the following relation between the unphysical component /i and the scalar density: 

h = -\d t W . (4.18) 

Comparing Eqs. ( |4.17| ) from the equal-time Feshbach-Villars approach with the corre- 
sponding Eqs. ( |4.7[ ) from the covariant approach using our energy- averaging procedure, one 
realizes that they have a very different structure. The difference is even more obvious when 
looking at the corresponding equations for the space-time densities which are obtained by 
integrating Eqs. (|4.17|) over momentum: 

dtp(x,t) + V-j(x,t) = 0, (4.19a) 
d t e(x,t) +V-P(x,t) - E(x,t)-j(x,t) - —^d t j -4sW(x,p,t) = 0, (4.19b) 

e(x)= J^{\( d t-l V *) +p2 + m2 ) W{x > P > t) - (419c) 
The first equation describes charge conservation and agrees with our Eq. fl4.9a|) . When 



comparing Eqs. ( 4.9b|) and ( |4.19b| ) one sees that the energy-momentum conservation law 



derived from the Feshbach-Villars approach features a strange additional term of unknown 
physical origin. Similarly, when compared to Eq. (|4.9cj) , Eq. (|4.19c| ) contains a strange 



factor | in front of the Laplace operator which breaks Lorentz invariance. We believe 
that these discrepancies indicate a technical problem with the Feshbach-Villars approach of 
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Ref. ||. Whether the problem resides in the derivation of the transport equations (|4. 14|) 
or in the postulated relations ( [4.15[ ) between the Feshbach-Villars spinor components /j 



and the physical densities could not be clarified. Only in the case of spatially constant 
electromagnetic fields (which to our knowledge is the only situation to which these equations 
have so far been applied [§|,|T2J) the two se ^ s °f equations Q4.9p and ( |4. 19| ) coincide. If we 
further neglect the magnetic field B, we derive two decoupled ordinary differential equations 
in time for W and p: 

d t p = , 

(d\ + 4(p 2 + m 2 )d t + AeE-p) W = . (4.20) 

These equations are identical in both approaches. The first equation expresses charge con- 
servation in a homogeneous electric field, and the second equation is the well-known equation 



of motion for the charge current j [11.121 which, according to the relation j = e , p — n, 

y/p 2 +m 2 

governs the time evolution of the particle density n by pair production in the electric field. 
This latter quantity thus comes out the same in both approaches. 

C. Spinor QED 

The case of spinor QED has been discussed in the context of the energy averaging method 
in Refs. |T^JT7|Jl~^| . The discussion presented in those papers has, however, focussed entirely 



on the equal-time kinetic equations for the lowest energy moment of the covariant Wigner 
function. Here we will reformulate the problem in terms of the equal-time hierarchy of 
energy moments as introduced in Sec. [TI| and also discuss the equations of motion for the 
higher order moments. 

Let us briefly review the relevant technical steps. We begin by performing a spinor 
decomposition || of the covariant Wigner function, separating in a second step explicitly 
the temporal and spatial parts of the covariant spinor components: 



W(x,p) = J 



F(x,p) + i"f S P(x,p) + lfl V^(x,p) + -f^ 5 A^(x,p) + K^ v {x,p) 
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(4.21a) 



j[loF (x,p) + 7 5 7oF 1 (x,p) + i^ 5 F 2 (x,p) + F 3 (x,p) 

-l5l-G (x,p) - 7-Gi(x,p) + i<y 5 j-G 2 (x,p) + 7 5 7o7-G 3 (x,p) . (4.21b) 



Inserting the decomposition (|4 . 2 1 a| ) into the covariant (VGE) equation of motion || for the 
Wigner function W(x,p), 

fr{%{x,p) + fafap)) - m) Wfop) = , (4.22) 

with Il^Xjp) and D^(x,p) from Eqs. ( |4.3| c,d) , and separating real and imaginary parts one 
arrives at two groups of coupled covariant kinetic equations of type ( |2.2|) , see Eqs. (74,75) 
of Ref. [ I2[ . Performing the energy average then leads to two groups of equal-time kinetic 



hierarchies (with 16 such hierarchies of equations in each group) for the energy moments of 
the 16 covariant spinor components. 

Since for spinor QED M = 0, minimal truncation of these 2 x 16 hierarchies according to 



Sec. [II B| results in one transport equation from each of the 16 transport hierarchies and no 
constraint equations. The minimal subgroup of equal-time kinetic equations thus consists 
only of the 16 transport equations for the 16 zeroth-order energy moments fi(x,p,t) and 
gi(x,p,t) of the covariant spinor components Fi(x,p) and Gi(x,p) {i = 0,1,2,3): 



dtfo + dg 1 = , 




(4.23a) 


dtfi + d-g + 2mf 2 = , 




(4.23b) 


d t f 2 + 27T-</3 - 2mA = , 




(4.23c) 


d t f 3 - 2n-g 2 = , 




(4.23d) 


dt9o + dfx - 2tt x g 1 = , 




(4.23e) 


d t gi + df - 2tt x g + 2mg 2 


= , 


(4.23f) 


d t g 2 + dxg 3 + 2tt/ 3 - 2mg 1 


= , 


(4.23g) 


d t g 3 - dx g 2 - 2nf 2 = . 




(4.23h) 



These 16 equations are identical with the BGR equations derived by Bialynicki-Birula, 
Gornicki and Rafelski @ by Wigner transforming the equations of motion for the equal- 
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time density matrix. They determine the dynamics of the zeroth-order energy moments. 

The three-dimensional dynamical operators occurring in these equations are identical with 

the ones arising in scalar QED and are given in Appendix 

The first-order moments fl(x,p) and gj(x,p) satisfy 16 transport equations derived (12 

from the first energy moment of the covariant transport equations ( |2.2a| ), 
1 



^d t fl + -j=d-g\ + D/ + I 9l = , 
-\ftfl + -^d-g l + A m /i + Dh + I go = , 
^m<kfl + ^-9l - ^fifl + Df 2 + 2Gg 3 = , 
~^d t fl - ^f-9\ + Dh - 2Gg 2 = , 

-J^al + -j=df\ + -^=tt x ^ + D 9o + Ifi ~ 2G x 9l = , 

1 * , 1 - , 2 , 2 , 

7/^*3i + ~^ d ^ ~ ~J% * X9 ° + 7f m ^2 + D 9i + If - 2G x c/o = 

1 * 1 - 2 t 2 , 

+ ^ d X 9 3 + - ~^ m 9i + D 9 2 + Ixg 3 + 2Gf 3 = 

1 - 1 - 2 " 

-^ d t9\ ~ xg\- -j=*fl + Dg 3 - I x g 2 - 2Gf 2 = , 



(4.24a) 
(4.24b) 
(4.24c) 
(4.24d) 
(4.24e) 
(4.24f) 
(4.24g) 
(4.24h) 



and 16 constraint equations derived [12] from the zeroth energy moment of the covariant 
constraint equation (|2.2b| ): 



^=/o = *-9i ~ TToio + mf 3 
~^fi = n-g - 7r /i > 



1 

i 

1 

71 

1 

7s 
i 

7^ 



s 



9l 



n 



d X flfx + 7T/i - TTofl'O + ™C/ 3 

dx g + irf - TT g 1 , 

^d/ 3 + 7T X £3 - 7r o 2 , 
-d/ 2 - 7T X flf 2 - 71003 + ^00 



(4.25a) 
(4.25b) 
(4.25c) 
(4.25d) 
(4.25e) 
(4.25f) 
(4.25g) 
(4.25h) 



28 



A discussion similar to that in scalar QED reveals |T8] that the transport equations (f4.24| ) 



are not independent of the BGR equations ( |4.23| ) and the contraint equations ( |4.25| ). For 



instance, using Eqs. ( |4.25| ), the transport equation for /q can be expressed as an operator 
combination of the transport equations for fo and g\\ 



7T 



(d t fo + d-flfi) - 7T- (d t g x + df - 2tt x g + 2mg 2 ) = . (4.26) 



Therefore, the first-order moments are fully determined in terms of the solutions of the BGR 
equations ( [4.23| ) for the zeroth-order moments by solving the constraint equations ( |4.25| ). 



Again, we have not been able to find a simple proof that the same is generally true for 
all higher order energy moments, and we stopped here. We do, however, believe that such 
a proof must exist |19[], and that therefore all higher order energy moments can be directly 
computed from the solutions of the BGR transport equations by solving the constraint 
hierarchy. 



It was shown in Ref. |L2| that in the classical limit the simple algebraic relation ( |1.4| 



changes the structure of the constraints ([4.25|) for the first-order moments and turns them 



into additional constraints for the zeroth-order moments (i.e. the equal-time Wigner func- 
tions). These extra constraints reduce the number of independent zeroth-order moments 
from 16 in the quantum case to 4 in the classical limit. They are thus extremely important. 
As a result the BGR equations reduce to two decoupled Vlasov-type transport equations for 
the charge and spin distribution functions. In the general quantum case there are no such 
extra constraints on the equal-time Wigner functions [|I8|| . One must solve all 16 coupled 
transport equations ( ^4.23| ), but these solutions then fully determine also all higher order 
moments. These higher order moments have important physical meaning: the first-order 
moment of Fo(x,p), for instance, describes the energy distribution in phase-space. With the 
help of the contraint ( |4.25a| ) it is given by 

e(x,p,t) = Try _£ 7o W(x,p) = ^-f^(x,p,t) 

= \/2~(mf 3 (x,p,t) -n f (x,p,t)+ir-g 1 (x,p,tfj . (4.27) 
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V. CONCLUSIONS 



We have presented a universal method for the construction of equal-time quantum trans- 
port theories from the covariant quantum field equations of motion. It is based on energy 
averaging the covariant kinetic equations for the covariant Wigner operator (which is the 
Wigner transform of the covariant, "two-time" density matrix) and its energy moments. 
This procedure yields a hierarchy of coupled transport and constraint equations for the en- 
ergy moments of the covariant Wigner function, the so-called equal-time Wigner functions. 
We showed how, in the mean-field approximation, this hierarchy can be truncated at a rather 
low level, requiring the solution of only a small number of equal-time transport equations, 
and how the higher order energy moments (higher order equal-time Wigner functions) can 
be constructed from these solutions via constraints. 

The major advantage of the equal-time formulation of (quantum) transport theory is 
that the resulting transport equations can be solved as initial value problems, with bound- 
ary values for the equal-time Wigner functions at t = — oo which can be calculated from 
the fields at t — — oo. This is not the case for the covariant transport equations and the 
covariant Wigner function. The present paper thus provides an essential step in the di- 
rection of practical computations of the dynamics of relativistic quantum field systems out 
of thermal equilibrium in the language of transport theory, i.e. in a phase-space oriented 
approach. The method presented here improves upon previous approaches by being much 
more systematic: we did not just focus on the lowest energy moments (which contain only 
a small fraction of the information contained in the covariant Wigner function), but we 
discussed and showed how to solve the complete hierarchy of moment equations. We had 
already before demonstrated for spinor QED that the non-covariant three-dimensional ap- 
proach (which starts directly from the equal-time density matrix) yields an incomplete set of 
equal-time transport equations. In this paper we also discussed the case of scalar QED and 
again discovered serious problems with the direct three-dimensional approach based on the 
Feshbach-Villars formulation. We conclude that the only safe way of deriving a correct and 
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complete set of equal-time quantum transport equations is by starting from the covariant 
formulation and taking energy moments of the covariant kinetic equations. The method 



can be generalized in a straightforward way to other types of interactions [|16|-|18|j , including 
non-Abelian gauge interactions fll9|| . 

The structure of the hierarchy of equal-time quantum kinetic equations depends on the 
structure of the covariant field equations from which one starts. For scalar or vector theories 
with second order time derivatives one has to solve a coupled set of three equations for the 
three lowest energy moments, two resulting from the equal-time transport hierarchy and 
one stemming from the constraint hierarchy. For spinor theories with only first order time 
derivatives in the field equations one ends up with only one equal-time transport equation 
for the lowest energy moment of each spinor component of the Wigner function. All higher 
order energy moments can be determined from the solutions of these equations by solving 
constraints. 

Important further simplifications occur in the classical limit h — > 0: then all higher order 
energy moments can be expressed algebraically in terms of the zeroth energy moment, and 
the number of equations is drastically reduced. For scalar theories one obtains just one 
Vlasov-type equation for the on-shell charge distribution function. For spinor theories one 
obtains two decoupled Vlasov-type equations (one scalar and one vector equation) for the 
on-shell charge and spin density distributions in phase-space. Again the only systematic 
way of deriving the constraints leading to these simplifications is by energy averaging the 
(classical limit of the) covariant transport equations. 

All results in this paper were derived in the mean field approximation, i.e. in the colli- 
sionless limit. It is generally known that including collision terms in the covariant transport 
equations leads to the appearance of non-localities in time ( "memory effects" ) in the equal- 



time transport equations ptifl . It is not inconceivable that these memory effects lead to 
serious complications for the truncation of the equal-time transport hierarchy. This is cer- 
tainly an interesting and difficult problem for future studies. 
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APPENDIX A: THE COEFFICIENTS <2? n 



^From the definition of the coefficients c^ n in ( 2.18Q and the orthogonality condition 
( p.5|) we see that 



c'S = % • (Al) 



Furthermore, it is easy to see from Q2.18 ) that 



c^ n = for n > i + m or j > i + m — n . (A2) 



Therefore, if the finite sum over m has the upper limit M, the sum over j in Eq. (p. 16 ) 
extends only over the range j < i + M. 

^From the recursion relation for the Legendre polynomials 

"PM = ^r^+iM + 2iVi Pi - l{uj) (A3) 

one easily derives the recursion relation 

m+l,n ^ "I - 1 m,n . ^ m,n ( A A\ 

^(2? + l)(2z + 3) v /(2z + l)(2z-l) 



for the coefficients c™ n . This allows to raise the first upper index m, starting from (AT). In 
particular we have 

- 10 - ' (A5a) 



(2i - l)(2i + l) 



£ w = 7 r4=== = , (A5b) 

(2i + l)^/(2z-l)(2z + 3) 
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which we will need in Appendix |B[ Similarly the second upper index can be raised by using 
the relations 



d u (P j+1 (u) - Pj-i(u)) = (2j + 1)P» , ^(±1) = (±iy (j > 0) , (A6) 
which, for j > 2, lead to 



This expression is useless for j = and j = 1 ; these cases can be treated by another recursion 
relation which can be obtained by using (|A6|) on the Legendre polynomial with the index i: 



m,n+l / 2i ~H 1 m,n+l I /<-> • \ 1 \ /r> • 7T m,n I m—l,n / 2z H 1 m —\ n \ / \ o\ 

For i < 2 this must be used together with 

cCC = v^TCt, (A9) 

which follows from P-i(x) = Pj_i(x). 

The above recursion relations can be initialized with the following nonvanishing coeffi- 
cients for m, n < 1: 



„oo _ „oo _ 11 
c oo ~~ c n ~~ c oo 



1, c i? = c ^ = i=, <# = ->/3, eg = -2. (A10) 



v/3 



APPENDIX B: THE OPERATORS g %j AND FOR SCALAR THEORIES 

The three-dimensional dynamical operators §ij(x,p) and fij(x,p) needed in Eqs. ( p.!4|) 
and ( |3.19| ) are obtained from the definition ( 2.17|) with the coefficients c^ n from Appendix 



Please note that for the upper limit in fl2.17]) for the sum over m is M = 1 while for 
it is M + 1 = 2. One finds: 

<?oo = Coo — , (Bla) 
001 = 7/^1° ' ( Blb ) 
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5io = V3 (-G01 + Ig 10 + 2G 12 } . (BID 



g u = G 00 ~ 2G n , (Bid) 
2 

=( 

15 



2 * 

5i2 = ^=^Ci , (Ble 



V^(3G 2-Gii-9Gi 3 ) , (Blf) 



520 

£ 21 = Vl5f-Goi + ^10 + 3612) , (Big) 



15 

522 = Goo - 3G U , (Blh) 
to = -Lg 10 , (Bli) 



35 

foo = Foo — Fn + -F 20 + 2F 2 2 , (Blj) 

/01 = (Ao - 2F 2 i) , (Blk) 

/02 = ^Ao , (Bll) 

Ao = V3 (-Ai + -Ao + 2Fi 2 - Ai - 6F 23 ) , (Blm) 

hi = ho ~ 2Ai + + 6A2 , (Bin) 

/ la = -|= (Ao - 3Ai) , (Bio) 



15 

As = j^Ao • (Blp) 



The dynamical operators G mn and F mn are given in equations ( 3.12 ) for scalar mean field 
interactions and in Appendix for scalar QED. Please note that in both cases A* = for 
k 7^ and thus some of the expressions for A above simplify. 
We also note that in 

2 i-l+m 

A-l,i+l = C T-l,i+lFran (B2) 

m=0 n=0 

the only term which doesn't violate the second inequality in ( |A2| ) is the one with m = 2, n = 
0: 



A-i,t+i - c i-i,i+iAo - , + ; Ao • (B3) 

(21 + 1)^/(21- l)(2i + 3) 
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In the last equality we used ( |A5b|) . From Eqs. ( [3.12b|) or ( |C4g|) we thus see that /j-i^+i is 
just a constant. Following the same reasoning one also obtains 



1 i+m 

9i-i,i = ^2 ^2 c™L n lti G mn = c i £ lji G f io 

m=0 n=0 



-G 
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(B4) 



\2i - l)(2i + 1) 

where according to Eqs. ( |3. 12a|) or (|C3d| ) G\q is proportional to the time derivative dt resp. 



dL. 



APPENDIX C: THE OPERATORS G mn AND F mn FOR SCALAR QED 



The elementary operators fl M and are extensions of the covariant momentum and 
the covariant derivative d^, respectively. Their double expansions in uo and d/duo are 



Il = AuO + 7T + 

A) = <*t 



Ad B d 2 



C d 3 



A duo A 2 duo 2 A 3 &</ 
D d Ed 2 F d 3 



+ ■ 



+ 



A duo A 2 duo 2 A 3 duo 3 

n = 7T - — — + — — + ■ ■ • , 



+ ■ 



D 



A duo A 2 duo 2 

+ XaJ + A2a7 + 



with the equal-time dynamical operators, 

i 

tto^, P) = ieh / ds s E(x + ish'V p ,t)-'Vp , 

"2 

7r(x, p) = p — ieh 2 ds s B(x + isKV p , i) x V p 

"2 

d t (x,p) = hd t + eh 2 ds E(x + ish'V p ,t)-'V p 



d(x, p) = H Vj + eh Ids B(x + ishV p , t) x V 5 



p ■ 



A(ar,p) = eh z i dss 2 —E(x + isftV p ,t)-V 3 



B(x,p) 
C(x,p) 



ieh 3 - l 



dt' 

d 



t s 3 —E(x + ishV p , t)-V p , 



eft 4 /■! , 4 (9 3 



ds s —^E(x + ishV v ,t)-V v 

1 at 3 



(Cla) 
(Clb) 
(Clc) 
(Cld) 



(C2a) 
(C2b) 
(C2c) 
(C2d) 
(C2e) 
(C2f) 
(C2g) 
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3 . d 



D(x,p) = ieh 2 i dss—E(x + ish^ p ,ty'Vp, (C2h) 
eh 3 rh d 2 

E{x,p) = — J^dss 2 —E(x + iS hV p ,t)-V p , (C2i) 

F{x,p) = — J_ i dss 3 —E{x + tshVp,t)-V p , (C2j) 

G(x,p) = eh 2 J 2 i dss 2 ^B(x + ishV p ,t) xV p + ieh J 2 ± ds sE(x + is%V p ,t) , (C2k) 
JET(a;, p) = -y- y x s 3 —B(x + isftV p , t) x V p 



2 

.1 

-eh 2 I 2 dss 2 ^-E(x + ishV p ,t), (C21) 

I(x,p) = ieh 2 f 2 ds s^-B(x + ishV p ,t) x V p — eh f 2 dsE(x + isKV p ,t) , (C2m) 
J-- ot J-~ 

eh 3 rh , d 2 



j ,9 



* , . en f2 , n C . _ , _ 

J(ar, p) = — y i rfs s Qp B { x + ishV p , t) x V p 



2 

9 



+rf / ' ds sj-E(x + zsftVp, £). (C2n) 

"a 



By substituting these expansions into the expressions of G and F in (|4.3| ) and using the 
definition of G mn ( [2.9Q , we list some low-order G mn and F mn which will be used in the 
derivation of the kinetic equations for low-order energy moments, 



Goo = ^odt + ft-d , (C3a) 
Goi = ~ (jt D + 7T-J - Adt + G-d) , (C3b) 
G 02 = (H-d + G-) - 7T J - M - AD - tt 0j e) , (C3c) 



A 2 

G 10 = Ad t , (C3d) 

Gn = -/) , (C3e) 

Gi 2 = -~ , (C3f) 

G 13 = ^ , (C3g) 



and 



F 00 = -id 2 - d 2 ) - (A + tt 2 - tt 2 ) + m 2 , (C4a) 
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*bi = J 1} - R D}) - 2n A - {tt, G} + 2s) , (C4b) 

F02 = ^ Q (-D 2 - / 2 + {d, J} - {1, E}) - (A 2 -G 2 - {tt, H} - 2n B) + 3c) , (C4c) 

Ao = -2Att , (C4d) 

F u = -2i, (C4e) 

A 2 = 2^, (C4f) 

Ao = -A 2 • (C4g) 
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